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Two-loop renormalization group equations in the standard 
model are re-calculated. A new coefficient is found in the (3- 
function of the quartic coupling and a class of gauge invariants 
are found to be absent in the /3-functions of hadronic Yukawa 
couplings. The two-loop /3-function of the Higgs mass param- 
eter is presented in complete form. 
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Analysis based upon renormalization group equations 
(RGEs) plays an important role in the study of physics of 
the standard model (SM) and beyond. Detailed analysis 
of RGEs confirmed the behavior of asymptotic freedom in 
QCD, and thus helped to establish a non-Abelian gauge 
theory for the strong interaction [1]. The runnings of 
coupling constants and mass parameters are crucial in 
global analysis of high precision electroweak experiments 
[2]. On the other hand, RGEs analysis extrapolated to 
extremely high energy provides a possible test for physics 
beyond the SM. For example, gauge couplings do not 
unify within the SM. This gives extra evidence against 
simple grand unification theories such as SU (5) with- 
out supersymmetry, in addition to the non-observation of 
proton decay. On the other hand, gauge couplings seem 
to unify at a scale ~ 2 x 10^^ GeV in the minimal su- 
persymmetric standard model, which can be interpreted 
as an indirect evidence for supersymmetry as well as uni- 
fication theories [3-5]. Comprehensive analysis can be 
found in [6]. 

Computations of RGEs in gauge theories have been 
performed for various models to different orders of per- 
turbation. Persistent efforts yielded recently a four-loop 
result of the /3-function of the strong coupling constant 
[7] . Two-loop RGEs of dimensionless couplings in a gen- 
eral gauge theory as well as the specific case of the SM 
had been calculated long ago in a series of classic pa- 
pers by Machacek and Vaughn [8-10]. By introducing 
a non-propagating gauge-singlet "dummy" scalar field, 
two-loop RGEs of dimensional couplings can be readily 
inferred from dimensionless results [11,12]. These were 
used to derive the RGEs of supersymmetric theories a 
decade later [11]. 

In this paper we re-calculate the two-loop RGEs in 
the SM, in a combination of using the general results of 
[8-10] and direct calculations from Feynman diagrams. 
A new coefficient is found in the /3-function of the quar- 
tic coupling and a class of gauge invariants are found to 
be absent in /3-functions of hadronic Yukawa couplings. 



We will also present the two-loop /3-function of the Higgs 
mass parameter in complete form, which provides a par- 
tial but useful check on the calculation of the quartic 
coupling. Whenever discrepancy with the literature ap- 
pears, we carefully inspect relevant Feynman diagrams 
to ensure consistency. 

To fix notations, we define Yukawa couplings and the 
Higgs potential in the SM to be 
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where three families of fermions are grouped together so 
F/,, Fjj, H are 3x3 complex matrices, and 0^ = iT2(t>* ■ 
For each coupling constant x in Eq.(l), we define a cor- 
responding /3-function 
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where (3^'^ , (3^^ denote the one-loop and two-loop contri- 
butions, respectively. We use dimensional rcgularization 
and the modified minimal subtraction scheme (MS) for 
renormalization. The expressions of the /Si^'^'s are quite 
standard which can be easily reproduced. The evaluation 
of /3i 's will be the object of this article. 

Following the conventions of [8-10], we define the fol- 
lowing combinations of Yukawa matrices for later conve- 
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The complex Higgs doublet has to be decomposed into 
real fields. Further complication arises since [8-10] as- 
sumed implicitly that the fermion fields are real, while 
usual Weyl fermions are complex. Caution should be 
taken when Yukawa couplings and gauge representation 
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matrices of fcrmions arc dealt with [12]. All issues taken 
into account, /3-functions in the SM can be obtained in a 
straightforward manner. The lengthy algebra is greatly 
simplified with the aid of the symbolic software FORM 
[13]. The /?- functions of the gauge coupling constants 
are readily reproduced and conform to those in the lit- 
erature [8]. Firstly, we present the /3-functions of the 
Yukawa couplings. To one loop, 
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to two loops, 

^~'(^H = ^(H+H)2 - H+HF+F,3 - If+FzjH+H 
+ ^(F+F^)^ +^2(5) (^^FlFn - ^H+h) 
-X4(5) + ^A2-6AH+H + ^y4(5) 

^9l - ^9l + I653) F^Fz, 



where SU{3) x SU{2) x U{1) gauge couphng constants 
gs, 52, and gi are normahzed based upon SU{5), so the 
standard electroweak gauge coupling constants g and g 
are related to those by, g^ = (/| and g ^ = 3/5(7^. The 
matrices F^, F^, and H do not have to be invertible. 
Their inverses should only be understood symbolically 
and need not be introduced in principle. Properly inter- 
preted, the /3-functions are equal to the Yukawa matri- 
ces themselves multiplied by the right-hand side of the 
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corresponding equations. For in [9], there was the 
term — 2AHF^F£); which is absent in Eq. (6). A close 
inspection indicates that this term arises only from the 
Feynman diagrams shown in Figs. 1(a) and (b). How- 
ever, these two diagrams cancel with each other, which 
can easily be verified by an elementary calculation. Sim- 
ilarly in Eq. (7) for /3^^, the term -2AF£,H+H is also 
absent, in contrast with [9]. The corresponding Feynman 
diagrams are shown in Fig. 1(c) and (d) and again they 
cancel with each other. The leptonic results are included 
here for completeness. 
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FIG. 1. Two-loop diagrams which affect hadronic Yukawa 
couplings, (a) and (b) ((c) and (d)) cancel with each other, 
thus resulting in a null contribution to /J^"* (0 



(2)-. 



187 
"80 
79 



■5i 



135 
l6 



52^ + 1651 F+Fo 



9i 



16 



92 



80 
29 1 
200 ^is''^ 



IQgl \ H+H 



9t 



27 Q 9 31 o o 



35 

T 



52 + ^9l9l 



9 



_ /404 



80 



4 

53 > 



(7) 



The /3-function of the quartic coupling can be obtained 
in a similar manner. The calculation is greatly simplified 
by the fact that there is only one independent quartic 
coupling in the SM. Calculation for models beyond the 
SM with numerous quartic couplings would be more in- 
volved, due to proliferation of combinatorics. The one- 
loop contribution to the /3-function of A is 
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to two loops. 
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Note in Eq. (10), the coefficient of the term 
ATr(H+HF+FD) is -42, instead of 6 as given by 
[10]. Wc note that terms proportional to \H{S) and 
ATr(H+HF^F£)) arise partly from scalar boson propa- 
gators. The relevant Feynman diagrams are shown in 
Figs. 2(a) and (b). There are also related two-loop 
proper scalar quartic vertex diagrams, which are generi- 
cally shown in Figs. 2(c) and (d). It turns out that Fig. 
2(c) does not contribute to the /3- function, so only Figs. 
2(a), 2(b), and 2(d) need to be evaluated. In addition 
to calculating Feynman diagrams directly, we compare 
the coefficient of \H{S) and that of ATr(H+HF+Fr,). 
By including all specific diagrams and carefully collect- 
ing all coefficients, we find that the ratio of the term 
proportional to the latter over the term proportional to 
the former is 42, instead of —6. This substantiates Eq. 
(10). On the other hand, coefficients of terms Xgfg2 and 
Xgf in Eq. (10) conform to those in [14], and in fact Eq. 
(10) agrees with the corresponding result in [14] if only 
the top-quark Yukawa coupling is retained. 

The /3-function of can be inferred from the results 
in [10] by introducing a non-propagating "dummy" real 
scalar field (j)d with no gauge interactions, and carefully 
computing the combinatorics associated with the sym- 
metry factor [11,12]. Specifically, the mass term can be 
re-written as 
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by decomposing the complex doublet (f) into four real 
scalars (jji, (i = 1,4). li (j)d is taken to have no other 
interactions, then the /3-functions of has the same 
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FIG. 2. (a) and (b): part of hadronic Yukawa coupling con- 
tribution to Higgs boson propagators, which in turn afltoct [3\ ; 
(c) and (d): relevant proper scalar quartic vertex diagrams. 



form as that of the newly introduced quartic coupling 
Xddij- To one-loop, the function is 
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and to two-loops 
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Note the terms Uggf and ngg2 in Eq. (12), which orig- 
inates from generic Feynman diagrams shown in Fig. 3. 
These diagrams are proportional to the fermion Dynkin 
indices in the pure fermionic loops, while Dynkin indices 
are always positive. So all fcrmions in the loop contribute 
additively, thus resulting in the Ug factor in the expres- 
sion. Eq. (12) is consistent with the newly revised result 
of [14]. 




FIG. 3. Part of SU{2) x U{1) contributions to the Higgs 
boson propagators. 

To a good approximation, we can express the SM ef- 
fective potential of the Higgs field 

K//(<A) = m\t)Z\t)cl,+cl, + ^-X{t)Z\t){cl,+cl,f, (13) 
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in term of the running coupling constants and the run- 
ning Higgs mass, 
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where the generic symbol x represents all dimensionless 
couplings, including the gauge couplings, the Yukawa 
couplings and the quartic scalar coupling; 7^ is the 
anomalous dimension of the Higgs field, to be found in 
[8,12]; t = log(/)/M, where M is an arbitrary mass scale, 
at which the initial values of the coupling constants and 
the mass are defined. The vacuum expectation value of 
the Higgs field is determined by minimization of Vg//. 
The physical mass of the Higgs boson is simply the sec- 
ond derivative of Vgff evaluated at the minimum [15]. 

In summary, we have re-calculated the RGEs in the 
SM. A new coefficient is found in the /3-function of the 
quartic coupling and a class of gauge invariants are found 
to be absent in /3- functions of hadronic Yukawa couplings. 
The /^-function of the Higgs mass parameter is also pre- 
sented in complete form. The changes in Yukawa cou- 
plings affect the running of the CKM matrix and the 
quark masses. The changes in and A will change 
the Higgs potential, which in turn affect the triviality 
and vacuum stability bound on the Higgs mass. Due 
to the dominance of one-loop results and relative bigger 
contributions from gauge couplings over those from the 
quartic and Yukawa couplings (with the exception of the 
top quark) , numerical changes are not expected to signif- 
icant. For A and Yukawa couplings related the &-quark, 
the changes are magnified by the factor mf/M^, but 
again suppressed by the factor m^/M^. All these shifts 
will be included in a future comprehensive analysis [12]. 
On the other hand, for a heavy fourth family of fermions, 
the changes would be sizable. 

Note added: Upon completion of this work, we were 
informed that part of the results presented here were also 
reached in [16]. 
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